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Abstract 

We study the variation of Lyapunov exponents of simple dynamical 
systems near attractor-widening and attractor-merging crises. The 
largest Lyapunov exponent has universal behaviour, showing abrupt 
variation as a function of the control parameter as the system passes 
through the crisis point, either in the value itself, in the case of the 
attractor-widening crisis, or in the slope, for attractor merging crises. 
The distribution of local Lyapunov exponents is very different for the 
two cases: the fluctuations remain constant through a merging crisis, 
but there is a dramatic increase in the fluctuations at a widening crisis. 



PACS numbers:05.45.+b,05.70.Fh 



I. INTRODUCTION 



In this paper, we study the behaviour of the Lyapunov exponent in systems 
where there are abrupt changes in the dynamics as a parameter is varied. 
Our interest is in exploring the typical dependence of the maximal Lyapunov 
exponent (MLE) on the control parameter so as to elucidate the signature of 
a transition in the nature of the dynamics. 

In the context of dynamical systems, abrupt changes in the phase space 
most commonly occur at the so-called crises |]J, which are caused by the 
collision of a chaotic attractor with the stable manifold of an unstable periodic 
orbit. The three major types of crises are distinguished by the nature of 
discontinuous change they induce in the chaotic attractor. At a boundary 
crisis, the chaotic attractor is suddenly destroyed and replaced by a chaotic 
transient as the parameter passes through its critical value. This occurs 
when the attractor collides with the stable manifold of an unstable periodic 
orbit that lies on its basin boundary. At an interior crisis, a sudden increase 
or decrease in the size of the attractor occurs when the stable manifold of 
an unstable periodic orbit lying within the basin of attraction of the chaotic 
attractor collides with it. At an attractor-merging crisis two or more chaotic 
attractors simultaneously collide with the stable manifold of an unstable 
periodic orbit lying on their common basin boundary which results in the 
merging of the attractors. 

The qualitative change in the dynamics at a crisis is reflected in the Lya- 
punov exponents. The case of a boundary crisis is not very interesting since 



1 



the Lyapunov exponent is either zero (if the transient leads to a periodic at- 
tractor) or takes a value characteristic of the chaotic attractor onto which the 
trajectory ultimately lands. The variation of the MLE at typical interior and 
merging crises is more dramatic, and in this paper we study these phenomena 
in a variety of simple model systems [§, |3|, |], |5|. In all the crises, there is 
a similarity in the dependence of the Lyapunov exponent on the control pa- 
rameter. Our major observation is that MLE has a characteristic behaviour 
which is, however, distinct for the attractor-widening and attractor-merging 
cases. For interior crises which terminate a periodic window the dependence 
of MLE on the control parameter is sigmoidal, with a large increase in fluctu- 
ations subsequent to the crisis. This abrupt increase in the MLE at interior 
crises has been observed before || |7|, §, |9], [K| in some studies oil — d and 
2 — d maps and flows. In contradistinction the MLE only has a "knee" at 
attractor-merging crises: after the crisis, the rate of change of the Lyapunov 
exponent decreases significantly. Again, in contrast to the attractor-widening 
case, there is no attendant increase in the fluctuations of the local Lyapunov 
exponents subsequent to the crisis. 

In the next section we describe the phenomenology of the behaviour of 
the MLE at crises in simple maps and other low dimensional dynamical sys- 
tems. Our results have relevance to studies of systems at a phase transitions, 
especially as a number of recent simulations of realistic systems have looked 
at the Lyapunov exponent, K-entropy and related quantities as a function of 
temperature or other control parameters [13, 21, 22|. These considerations 
are discussed in relation to the present work in the concluding Section III. 
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II. The Lyapunov exponent at Crisis 

The Lyapunov exponent, which is used to characterise the degree of chaoticity 
of a dynamical system gives the average rate of exponential divergence of 
two nearby trajectories |20]. In a n— dimensional dynamical system there 
are n Lyapunov exponents and the system is chaotic if at least one of them 
is positive while for regular dynamics all Lyapunov exponents are zero or 
negative. We focus on the largest of these, which is most simply defined as 

X m = lim X m (t) = lim ]- lira log (1) 

t->°c t^oo t d(0)^0 d(0) 

where d(0) is the initial separation between two trajectories, and d(t) is their 
separation after time t. A number of methods have been proposed in the 
literature to compute one or more of the Lyapunov exponents JT^, |T3| , [TJ 



Here we use the tangent space method [JT^] which is sufficient since we are 
interested primarily in the largest Lyapunov exponent. 

To study transient objects like repellers or semi-attractors one can also 
analyse the finite-time exponents X m (t) which are also defined in Eq. [l|. The 
instability fluctuations on an attractor can also be studied by dividing a 
long ergodic trajectory in segments of size t and calculating the Lyapunov 
exponent X m (t) for each of these. The probability density P(X m (t)) of the 
distribution of local Lyapunov exponents has the scaling form, for t — > oo, 

P(X m (t)) ~ exp(-^(A m (*))), (2) 

where ip(X m (t)) is a concave function with its minimum equal to zero at X m = 
A m (oo). ]TR ITM]. In a highly mixing system, the time correlations of X m (t) 



can be ignored and then by the central limit theorem P(X m (t)) is a gaussian 
and if>(X m ) is parabolic. However at crises points the gaussian distribution 
breaks down and ip(X m (t)) develops a cusp at its minimum jlR [T7| ]. 

A. Interior or widening crises 

We first consider the logistic map 

Xn+i = rx n (l — x n ) (3) 

It is well-known 0] that as the parameter r is increased the logistic map 
undergoes period-doubling cascade terminating at the accumulation point 
Too ~ 3.5699 .... Beyond that the dynamics is mainly chaotic, punctuated 
at various intervals by periodic windows of arbitrarily high period. An odd 
period-n window is created at a saddle-node bifurcation together with an un- 
stable period-n orbit. As r is increased there is a sequence of period-doubling 
bifurcations creating periodic attractors of period 2n, 2 2 n, 2 3 n, .... Beyond 
the accumulation point of the period-doubling bifurcations the attractor is 
made up of n distinct pieces. The trajectory hops among these pieces in a 
regular manner but the distribution of points within each piece is random 
on the so-called semi-periodic attractor |Tj|. At the right end of the win- 
dow there is an interior crisis when each piece of the semi-periodic attractor 
meets a point of the unstable period-n orbit which was created in the saddle- 
node bifurcation, leading to an abrupt increase in the accessible phase space 
volume Q. 

Shown in Fig. 1 is the variation of A m with r — r™ near the n = 3, 5 and 
7-band crises, which occur at parameter r = r" respectively. We observe 
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that for all interior crises the MLE vs. (r — r") curve is sigmoidal. The 
standard deviation in the local Lyapunov exponent (calculated from iV = 10 4 
with 50 different initial conditions) increases dramatically at the crisis. It 
is easy to see why the fluctuation in the local exponents should increase 
abruptly at the interior crises: the attractor gains large volume which may 
have entirely different stability properties [K| [0| . The MLE increases at the 
crisis because the attractor engulfs the coexisting repeller. This repeller is 
the remnant of the chaotic attractor which had ceased to exist at the saddle- 
node bifurcation. Computation of finite-time Lyapunov exponent near crisis 
shows that the repeller has larger finite-time Lyapunov exponent than the 
semi-periodic attractor 0. The spectrum of local Lyapunov exponents of 
the post-critical attractor just before and after the 3-band crisis is shown in 
Fig. 2(a). The linear segments indicating non-hyperbolicity at the crisis are 
present on both sides of r c . After the crisis a kink at large A is visible which 
corresponds to the distribution on the repeller fT5|j . 

Pompe and Leven [[J argued that the increase in MLE is proportional to 
the probability density on the repeller, and model it by the power-law 

Pr ~ (r - r c y (4) 

We confirm that the increase in MLE is proportional to the probability 
density on the repeller. A power law fit to MLE data gives the exponent 
\x = 0.51 ± 0.04 for 3-band crisis and \i = 0.52 ± 0.04 for 5 and 7 band crises. 



Indeed, Grebogi et al. [pij obtained an approximate scaling near 
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^ (r _ r 3)l/^ (ln(r _ r 3 )); (5) 



where g is a periodic function. It therefore appears that this scaling relation 
is valid at all other band-crises as well. 

Other 1 — d maps also show the same phenomenology. We have studied 
a map originally introduced by Kariotis, Suhl and Eckmann|| in order to 
mimic the dynamical behaviour of intramolecular processes and isomeriza- 
tion. This map is given as 

x n+1 = rx n (u 3 -2ujx 2 n + x*), (6) 

We fixw = 0.8 and consider r as the control parameter. The above map shows 
both attractor-merging (discussed below) and attractor-widening crises. At 
the 5-band crisis at r\ = 5.2505109. . .., as shown in Fig. 1, clearly shows 
that the behaviour of the Lyapunov exponent is essentially identical to that 
observed for the logistic map. 

Higher dimensional systems also show the same behaviour: for example, 
the well known Henon map, 

y n+1 = bx n (7) 



which has a well-characterized, complex structure of bifurcations and crises pO 
Fixing 6=0.3 and varying r, the 7-band crisis occurs at rj=l. 2716856. Again 
(cf. Fig. 1) it is seen that the maximal Lyapunov exponent shows the by 
now familiar characteristic sigmoidal behaviour as the function of (r — rj). 
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B. Attractor-merging crises 

Attractor-merging crises typically occur because of some symmetry in the 
underlying dynamics, for example in the logistic map, beyond the accumula- 
tion point of the period-doubling cascade there is a successive merging of 
chaotic bands. Thus, for r < r m the chaotic attractor consists of 2 m chaotic 
bands. If we take a point in any one of these 2 m bands the trajectory will 
come to that band after 2 m iterations. So that the band can be considered as 
an attractor for 2 m -times iterated map. So the band- merging phenomenon 
can be regarded as an attractor-merging crisis KJ for the 2 m -times iterated 
map. 

We show results for the m — 1 merging at r m = 3.678486 . . . and the 
m = 2 merging at r m = 3.5925663 ... in Fig. 3. For band-merging crises 
the MLE vs. (r — r m ) curve has a sharp knee precisely at r = r m , i.e., 
the derivative of MLE is discontinuous at r c . However the knee angle is 
not same for all band-mergings. The local Lyapunov exponents are more 
uniform here, and we do not see any significant change in the fluctuation 
properties after the crisis. This is not difficult to understand as the co- 
merging attractors are symmetry-related. Similar behaviour was observed for 
attractor-mergings in the other maps studied |3|, |J (Fig. 3). The spectrum of 
local Lyapunov exponents near the attractor-merging crisis (Fig 2(b)) shows 
the linear segments indicative of the non-hyperbolicity at the crisis [[HJ . 

We also consider a merging-crisis in the forced Duffing equation |lj 

d 2 x/dt 2 + v dx/dt + ax 3 — [3x = r sin ut (8) 
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We take v = l,a = 100, (3 = 10, u = 3.5 and study Eq. (||) near its crisis 
value r m fh 0.853. Below r m there are two chaotic attractors one confined to 
the well in x > and another confined to the well in x < 0. These attractors 
merge at r = r m , where the MLE near r = r m again shows a knee (Fig. 3). 

III. Conclusions 

In this work we have studied the variation of the largest Lyapunov exponent 
MLE near interior and attractor-merging crisis for some simple and well- 
known systems. 

We observe that around an interior crisis, the MLE vs. the control 
parameter curve has a sigmoidal dependence, with the fluctuations increasing 
dramatically at the crisis. For the 3-band crisis in the logistic map (Fig. 1), 
for example, the average fluctuations just before the crisis is 4.2 x 10 -4 , while 
after the crisis it is 1.2 x 10 -3 . On the other hand, a knee-shaped curve is 
observed for attractor-merging crises with no increase in fluctuations beyond 
crisis. 

In recent work Fan and Chayflll]] have studied the Lyapunov exponents of 
Rose-Hindmarsh system (consisting of three coupled differential equations), 
and report that Lyapunov exponents are not good indicators of an interior 
crisis. They prefer the use of topological entropy which showed an abrupt 
increase at the crisis. However, in contrast to the present systems where the 
interior crises terminate a periodic window the interior crisis they studied 
was caused by collision of two period-adding bifurcation processes travelling 
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in opposite directions in the parameter space. This may be one reason why 
they did not observe increase in the MLE at the crisis. 

The observations made above are of relevance to recent simulation studies 
of systems undergoing a change in bulk phase. In recent years it has become 
possible to study the detailed dynamics of mesoscopic systems undergoing 
phase change, and a number of studies |14|, 21], ^] have therefore focussed 



on the relation between phase transitions and the Lyapunov exponents that 
characterise the dynamics. For example, in a study of a large number of cou- 
pled planar rotors, Butera and Caravati[{21]j found a discontinuity in the slope 
of the maximal Lyapunov exponent (MLE) at the precise temperature of the 
Kosterlitz-Thouless transition. More recently, it has been seen in molecu- 
lar dynamics simulations of small Lennard- Jones clusters p2| that the largest 
Lyapunov exponent also increases dramatically as the system makes a tran- 
sition from a solidlike to a liquidlike phase. This abrupt change characterizes 
the phase transition in a manner exactly analogous to the Lindemann cri- 
terion, and indeed offers an alternative connection between the phase space 
dynamics and the phase change. 

This abrupt increase in the MLE corresponds to an increase in the avail- 
able phase space volume and consequently in the local rate of divergence of 



trajectories[12|. Berry and co-workers [ 14] have looked at a variety of dynam- 
ical indicators, including the Kolmogorov-Sinai (KS) entropy, i.e., the sum 
of all positive Lyapunov exponents. This quantity increases smoothly with 
temperature or energy as the phase changes although information can be 



obtained regarding underlying potential-energy surface [14 
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The present study further underscores the utility of the MLE as an indi- 
cator of phase transformation since the band-crises between different chaotic 
phases are, in a sense, dynamical system analogues of phase transitions. Prior 
to the crisis there is long-time correlated noisy periodicity while after the cri- 
sis, dynamics lacks long-time correlation since motion is chaotic within a 
single-band attractor. 
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Figure Captions 

Fig. 1 Variation of the maximal Lyapunov exponent , X m , vs. (r — r c ) around 
interior crises. The error bars show the magnitude of the fluctuation 
in the local Lyapunov exponent. Beyond the crisis, solid line through 
the data is the power-law r — r c ) M , with \x = 0.5 for a), b) and c), 
and /j, = 0.44 for d) and 0.37 for e). 3-band crisis in logistic map at 
= 3.8568007. (b) 5-band crisis in the logistic map at r\ = 3.7447104. 
(c) 7-band crisis in the logistic map at r 7 c = 3.70279404. (d) 5-band 
crisis in the Kariotis-Suhl-Eckmann map at r\ = 5.2505109. (e) 7-band 
crisis in the henon map at rj = 1.2716856, b is fixed at 0.3. 

Fig. 2 Spectrum of local Lyapunov exponents, ip(X(n)), just before and after 
a crisis, for logistic map (a) the 3-band crisis- r = r% ± 5r . (b)the 
m = 2 band-merging- r = r m ± Sr. Here n = 60 and the number of 
iterations of the map is 3 x 10 6 . 5r = 1 x 10~ 6 . Crosses refer to data 
before crisis and circles to data after crisis. 

Fig. 3 Variation of the maximal Lyapunov exponent, X m , vs. (r — r c ) around 
attractor-merging crises. As in Fig. 1 the error bars show the magni- 
tude of the fluctuation in the local Lyapunov exponent. The straight 
lines through the data are least-square fit constrained to pass through 
the critical point, (a) m — 2 band- merging in the logistic map at 
r = 3.678486. (b) m — 1 band-merging in the logistic map at r = 
3.59256296. (c) Attractor-merging in Kariotis-Suhl-Eckmann map at 
r = 5.74009035. (d) Henon map- merging crisis at r = 1.084404, b is 
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again fixed at 0.3. (e) Well- merging in the forced Duffing equation at 
r « 0.853. 
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